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ABSTRACT 

A variety of physical heating mechanisms are combined with radiative cooHng to explore, via one dimen- 
sional hydrodynamic simulations, the expected thermal properties of the intracluster medium (ICM) in the 
context of the cooling flow problem. Energy injection from type la supernovae, thermal conduction, and 
dynamical friction (DP) from orbiting satellite galaxies are considered as potential heating mechanisms of the 
central regions of the ICM, both separately and in conjunction. The novel feature of this work is the exploration 
of a wide range of efficiencies of each heating process. While DP and conduction can provide a substantial 
amount of energy, neither mechanism operating alone can produce nor maintain an ICM in thermal balance over 
cosmological timescales, in stark contrast with observations. Por simulated clusters with initially isothermal 
temperature profiles, both mechanisms acting in combination result in long-term thermal balance for a range of 
ICM temperatures and for central electron densities less than ~ 0.02 cm"^; at greater densities catastrophic 
cooling invariably occurs. Furthermore, these heating mechanisms can neither produce nor maintain clusters 
with a declining temperature profile in the central regions, implying that the observed "cooling-core" clusters, 
which have such declining temperature profiles, cannot be maintained with these mechanisms alone. Super- 
novae heating also fails to maintain clusters in thermal balance for cosmological timescales since such heating 
is largely unresponsive to the properties of the ICM. Thus, while there appears to be an abundant supply of 
energy capable of heating the ICM in clusters, it is extremely difficult for the energy deposition to occur in 
such a way that the ICM remains in thermal balance over cosmological time-scales. Por intracluster media 
that are not in thermal balance, the addition of a small amount of relativistic pressure (provided by e.g. cosmic 
rays) forestalls neither catastrophic heating nor cooling. These conclusions are driven largely by the fact that 1) 
DP heating scales approximately as the gas density, while cooling scales as gas density squared, and thus DP 
heating cannot generically balance cooling without fine-tuning; 2) conduction acts to erase temperature gra- 
dients, while most observed clusters in fact show strong gradients in the inner regions. These results strongly 
suggest that a more dynamic heating process such as feedback from a central black hole is required to generate 
the properties of observed intracluster media. 

Subject headings: cooling flows — galaxies: clusters — hydrodynamics — conduction 



L INTRODUCTION 

Groups and clusters of galaxies are filled with hot plasma 
in at least approximate pressure equilibrium with the grav- 
itational potential of their dark matter halo. Observations 
have demonstrated for decades that the cooling time of cen- 
tral regions of this intracluster medium (ICM) in most (> 
70%) clusters is shorter than a Hubble time (e.g. Edge et al] 
T99alPeres et al.ll998l:[Sanderson et al.l2006tlvShlinin et all 



which carries heat from the abundant thermal reservoir 
of the cluster gas to the cooling inner parts, gravitational 
heating, including the heating due to the or bital motions 
of galaxies, i.e . dynamical friction (DP; e.g. lMillerlll986l: 



200637 indeed, it is often as short as ~ 0.1 - 1 Gyr. In the 
absence of heating, the ICM will thus cool and flow into 
the central galaxy at the prod igious rate of hundreds of so- 
lar m a sses per year (ICowie & Binnev 1977; Pabian & NulsenI 
119771: iMathews & Bregmanlll978i) ; see fabian, a994 for a 
review. 

Since this (catastrophically) cooling gas has never been ob- 
served (Peterson et al. 2001, 2003; Tamuraetal. 2001), it is 
now generally supposed that there exists one or more heating 
mechanisms that maintains the ICM in overall thermal bal- 
ance. The tension created by the fact that the ICM is not 
observed to be significantly cooling, despite the short cool- 
ing times near the cluster center, has become known as the 
cooling flow problem. 

Many possible heating mechanisms have been i nvesti- 
gated , including therma l con duction (e.g. [ Binney & Cowie| 



Just et al.lfT9 9^. 'Pabian"200 3t lEl-Zant et al.ll2()04btTKim'etan 
2005t iDek el & Birnboim n2007h. turbulent mixing (e.g 
^eiss & Just 1996; Kim & Nai-avan| |2003bl IVoigt & FabianI 
2004; Dennis & Chan dran 2005), energy injection from 
activ e galactic nuclei (AGN) via bubbles, sound waves, 
etc. (iBinnev & Taboi^' 19951; 'Nusser et al."2006'; 'Binney et alj 
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20071: Ruszkowski et al. 2004a b; Pujita & Suzuki 2005| 
Mathew s et al. 2006; Ciotti & Ostri kert2001tlVoit & Donahui 

2005t Ciotti & Ostriker 2007), and combinations of 

conduction and AG N (Ruszkowski & Begelman !2002l: 
iBrighenti & Mafliewsl 12003; .Puiita & Suzuki ,20051 
conduction and cosmic -rays (CRs; iGuo & OhI l2007h . 



or 



1981; Tucker & Rosner 1983; Voigt et al. 2002; Pabian et al.' 
Tmi; Zakamska & Nai-avan 2003; Kim & Naravan 2003a; 
Voig\ & FabianI 120041: iDolag et al.1 l2004t iPope et all l2006l) 



AGN and preheating ("M cCarfliy et al.l |2007a'). One of the 
most comprehensive studies to date was undertaken by 
iBrighenti & Mathewsl (l2002h who included many of the 
heating mechanisms mentioned above and found that no 
combination generically reproduced observations. A strin- 
gent constraint on the potential heating mechanism is that it 
must not only supply of order the necessary energy but must 
also supply it in a way that approximately balances cooling 
(locally) so that the ICM does not either heat up or cool on 
cosmological time-scales. 
Many of these mechanisms have been idealized, and it is 
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unclear how effective they would be in observed clusters. 
For example, though difficult to constrain observationally, 
it has been suggested that the conductivity in clusters may 
be much smaller than is required to st em the cooling flows 
dMarkevitch et al.'2003; Xiang et al. 2007). In addition, while 
DF heating may be attractive, since there are a plethora of 
satellite galaxies in clusters, it has yet to be identified as 
impor tant in hydrodynamical simulations (iFaltenbacher et all 
l2005h. although t his may be due to insufficient resolution 
(|Naab et al. 2007). Finally, bubbles from activity in cen- 
tral galaxies m ay be important (e.g. I Churazov et al.l 120011 
iRuszkowski & B egelman 2002; Briiggen & Kaiser 2002), but 
the mechanism by which they transfer their energy to the clus- 
ter gas remains obscure. 

Even if one or more of these mechanisms could supply the 
requisite energy (i.e. result in thermal balance), the mecha- 
nism must, in addi t ion, not allo w the gas to be thermally un- 
stable ( lFieldlll965l: lBalbus|[T986) . In what follows we make 
extensive use of this distinction between thermal balance, 
where the net cooling of a given parcel of gas is zero, and 
thermal stability, which pertains to the ability of a gas to re- 
main in thermal balance in the presence of isobaric pertur- 
bations. While attention in the literature has focused on in- 
vestigating the stability of an ICM in thermal balance (e.g. 
iKim & Naravan 2003b), herein we address the more funda- 
mental challenge of maintaining a cluster in thermal balance 
with one or more of the heating mechanisms described above. 

The physical mechanism causing thermal instability is easy 
to understand. Small regions which are slightly over-dense 
will radiate more than their surroundings (since cooling per 
unit volume scales as p^) and, isobarically contracting, find a 
new equilibrium, which is at a higher density and lower tem- 
perature. For the temperature domain in question, this leads 
to further cooling and a thermal runaway results. If a heating 
process such as conduction or DF heating exists, there will 
be a stabilizing influence. But if the process scales as p (per 
unit volume) then the equilibrium will not be stable, so high 
temperature under-dense regions will heat exponentially and 
lower temperature over-dense regions will still cool exponen- 
tially. However, a non-thermal component (such as cosmic 
rays or tan gled magnetic fields ) can in principle suppress this 
instability (ICenll20()5l; IGuo & Oh 2007, see also the appendix 
herein) because the non-thermal pressure partially decouples 
the hydrodynamic balance from the thermodynamics of the 
cluster gas. In other words, an ICM with a non-thermal com- 
ponent that radiatively cools will lose less total pressure sup- 
port than an ICM with no non-thermal component. An ICM 
with a non-thermal component will thus have to contract less 
to compensate for the lost thermal pressure support. 

Observed clusters can be rather cleanly divided into two 
categories based on the properties of their ICM. "Cooling 
core" (CC) clusters are those whic h have steep tempera- 
ture drops in their centr al regions ([Sanderson et al.l |2006|) 
and metallicity gradients (IDe Grandi & Molendil200l1) : many 
CC clusters have identified radio emission at their centers 
jBest et al.l l2005h and show si gns of AGN activ ity, includ- 
ing observed sound waves (e.g. lFabian et al.l l2006) apparently 
emanating from their center, and bubble s at an average dis - 
tance of 20 kpc from the cluster center (iBirzan et al.ll2004) . 
In contrast, non-c ooling core (NCC) clu sters are approxi- 
mately isothermal (ISanderson et alJl2006|) within ^ 100 kpc 
and s how little metallicity gradient (iDe Grandi & Molendil 
1200 lb . CC clusters generally have central cooUng times of 
0.1-1.0 Gyr while the cooling times in NCC clusters are gen- 



erally somewhat higher at ~ 1 Gyr ( ISanderson et al.ll2006h . 
The present work investigates the expected thermal properties 
of both CC and NCC clusters. 

The aim of the present study is not to find one or more heat- 
ing mechanisms that can offset radiative cooling but rather to 
understand the general conditions that lead to thermal balance 
in the ICM for a variety of possible mechanisms. This re- 
quires both a heating mechanism (or mechanisms) that can en- 
ergetically offset radiation, and also maintain long term ther- 
mal balance. 

The rest of this paper proceeds as follows. In f|2]the meth- 
ods are discussed, including the implementation of cooling, 
supernovae heating, conduction, DF heating, and relativis- 
tic pressure, the initial equilibria, and the numerical setup. 
Sj3] contains the results of a series of numerical experiments 
where the thermal balance of the ICM is investigated. A dis- 
cussion and summary of these results can be found in §f|4land 
|5] respectively. Throughout we assume h = 0.7 where h is the 
Hubble parameter in units of 100 km s"' Mpc"' . 

2. METHODS & PHYSICAL PROCESSES 

This section reviews the relevant fluid equations, cooling 
function, and heating mechanisms that will be explored in 
depth in the following sections. This section also discusses 
our implementation of a relativistic fluid component, initial 
equilibria, and the numerical setup. 

2. 1 . General Equations 
The hydrodynamical equations are: 



dt 



1 

Pg 



-VP,o, + g, 



deg 
~di 



-(e„ + R)V-v+r-A, 



(1) 



(2) 



(3) 



where pg, eg, Pg, ftot and v are the gas density, internal energy 
density, gas pressure, total pressure, and velocity, g is the to- 
tal gravitational acceleration, and T and A are the heating and 
cooling functions per unit volume. We have explicitly distin- 
guished between the total and gas pressures because below we 
will allow for the addition of a relativistic fluid that can pro- 
vide additional pressure support. The equation of state for the 
gas is: 

(4) 

where 7g = 5/3 is the ratio of specific heats. 

The gravitational acceleration is the combination of a (pas- 
sive) dark matter halo, a central cD galaxy, and the self- 
gravity of the gas: 



g(r) = g(r)r,M+8(r)cD- 
where is 



GMg,,« r) 



g(r) 



2GMn 



DM ■ 



ln(l + x) 



1 



(5) 



(6) 



x{l+x) 

and is derived from the NFW dNavarro et alj|1997h density 
profile: 

_ Mo/2n 
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where G is Newton's constant, x = r/r^, Mq is the normal- 
ization and Tj is the scale radius, i.e. the radius at which the 
density profile scales as r~^. For reference, the mass within 
2rj is equal to 0.86M() for the above density pro file. We fix 

= 460 kpc for comparison to iKim et al.l (l2005|). This value 
for is bracketed by the observational range (IVikhlinin et al.l 
|^06b). M() is allowed to vary as discussed in mA\ 

The central cD galaxy mass density profile is taken to be a 
King profile: 

PcD(r) = TT—T-, oTU? ' (°) 



[l + (r/reD)2]3/2' 



where 



PcD.O ■ 



(9) 



In the above equations, a^u is the central ID velocity disper- 
sion of the central galaxy and r^D is the core radius for the 
King profile; these parameters are taken to be 200 kms~' and 
2.83 kpc, respective ly, from a fit to the cD galaxy NGC 6166 
dKelson et al.ll2002h which is representative of cD stellar den- 
sity profiles. Note that the King profile is for all practical pur- 
poses quite similar to the more conventional de Vaucouleurs 
profile, but is more analytically tractable. The gravitational 
acceleration associated with this density distribution is 



CT^p /l + r'(l + r'2) 



./2\-0.5 



ln(r' + Vl + r'2) 



(10) 



gcD(r) = - 1 , 

'"cD V /(r' + Vl + r'2) 

where r' = r/r^D- 

Mass that flows through the inner boundary (1 kpc; see be- 
low), and is hence no longer followed directly in the simula- 
tion, is added to the cD galaxy. Our results are unchanged if 
the mass is added to a central black hole instead. 

The gas is assumed to be ideal: 

where ^ is the mean molecular weight, is the atomic mass 
unit, ks is Boltzmann's constant, rig is the electron number 
density and T is the gas temperature. Throughout, the gas 
metallicity is taken to be (1 /3)Zq (except in i)2. 41 where we 
discuss the sensitivity of our results to different metallicities). 

2.2. Heating & Cooling 

This section describes the various heating and cooling 
mechanisms that are potentially relevant for the thermody- 
namics of the ICM. Each mechanism includes an adjustable 
free parameter that is meant to encapsulate both our ignorance 
regarding the applicability of the mechanism to real clusters 
and the uncertain values of the particular parameters relevant 
for each mechanism; a summary of these parameters is pro- 
vided in Table [1] 

2.2.1. Radiation 

A simplified cooling function is adopted that broadly cap- 
tures the metal licity dependent features in the detailed cool- 
ing function of I Sutherland & Dopital (Il993h via the following 
function: 



A = 2.1 X 10 



-''nlT 



1/2 



1 + 



1.3 X 10^C(Z)y '^"' 



(12) 



in units of ergs" 'cm"-'. The variable ( is a simple function 
of the gas metallicity, Z. This metallicity-dependent cool- 
ing function is shown in Figure [T] (solid lines) along with the 
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Fig. 1 . — Cooling function used in the simulations (solid line s) com pared 
to the more accurate cooling functions of Sutherland & Dopita' U993I) . The 
solid lines are, from bottom to top, for primordial metallicity, [Fe/H] = -1.5 
and [Fe/H] = 0.0. Note that even for low metallicity the cooling function 
deviates from pure free-free cooUng (where A oc r'''-) at T > lO' K due to 
recombination cooling. 



detailed cooling functions from I Sutherland & Dopital (Il993h . 
The cooling function is truncated atT = 0.01 keV because our 
analytic approximation becomes inaccurate at lower tempera- 
ture, and because more complex physical phenomena not in- 
cluded herein (such as star formation) become relevant at such 
temperatures. Most of our results focus on a single metallicity 



of Z = ^ Z0 ; in !12.4l we briefly discuss how the results change 
(11) when different metallicities are adopted. 



2.2.2. Supernovae 

The stellar populations of central galaxies are extremely 
old, w ith typical formation epochs at z > 3 (iThomas et al.l 
120051) ; there are thus no type II supernovae events in these old 
systems. However, type la supernovae events are common, 
even for the old stellar populations comprising cD galaxies. 
At low redshift, roughly one type la super novae event oc- 
curs e very 100 years per 1O'^M0 of old stars dMannucci et al.l 
I2OO5I) . Since each supernova releases ^10 ergs, which we 
assume is transfered entirely to the ICM, this corresponds to a 
time-averaged energy injection rate of 10"^^ erg yr"' per galaxy 
of mass IO'^Mq. The type la rate for a co-evolving stellar 
population may be a decreasing function of time, indicating 
that the rate at earlier epochs would be higher than the value 
we adopt here (Greggio 2 0051 ; iMannucci et all2006) . For our 
purposes we are primarily interested in type la rates at > 1 
Gyr after the formation epoch of the stars, pre cisely where the 
time- dependence of the rates are least certain (IMannucci et alj 
I2OO6I) . However, even a modest increase of the rate with 
redshift would not qualitatively change our conclusions. We 
therefore do not include any time-dependence in the type la 
rate in what follows. 

Since there are hundreds of galaxies within clus- 
ters, supernovae-related energy injection should be treated 
throughout the cluster However, since radiation is capable 
of substantially cooling only the inner cluster region (because 
cooling scales as at fixed temperature), and since the ther- 
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mal reservoir in the outer regions is large, we only distribute 
supernovae within the central galaxy for simplicity. The su- 
pernovae energy injection rate is distributed with the same 
space density as the cD galaxy (cf. Equation[8]l. 

Supernovae energy injection is non-negligible in the cen- 
tral regions at the beginning of our simulations. For a initially 
isothermal cluster with central electron density 0.020 cm"^ 
and 7} = 6 keV, the volume-averaged supernovae heating ini- 
tially dominates cooling within only 2 kpc, while for a central 
density of 0.005 cm"^ supernovae heating dominates within 
~ 10 kpc. However, since the injection rate is independent of 
the properties of the ICM, it has little impact on the long-term 
thermal balance of the ICM. In particular, while for very low 
density intracluster media supernovae heating is comparable 
to radiative cooling on small scales, it is completely ineffec- 
tive at larger scales (i.e. several tens of kpc) because there 
are so few stars, and hence so few supernovae, at these larger 
scales. 

We have run many simulations for a variety of initial con- 
ditions and additional heating sources (discussed below) with 
supernovae feedback and indeed find that while at the central 
regions (~ 1 - 10 kpc) it can be important, it has little effect on 
the long-term thermal balance of th e ICM. These conc l usions 
are qualitatively similar to those of iKravtsov & YepesI (l2000h 
who used the observed metallicity of the ICM to demon- 
strate that supernovae heating is insufficient to offset radiative 
losses. In order to simplify the discussion in the following 
sections, we neglect this energy injection process for the re- 
mainder of this work. 

2.2.3. Conduction 

The heat flux due to electron conduction may be described 
by: 

Tcond = V-(AtVr), (13) 

where k is the conductivity. For a fully ionized plasma, the 
conductivity is: 



K=fKSp=f 



1.84 X 10-5 7^/2 
InAr 



erg s ' K ^ cm 



(14) 



where Ksp is the classicalfS pitzerl (1 1 962h conductivity, InAc ^ 
37 is the Coulomb logarithm, and / is a free parameter de- 
scribing the suppression of conductivity relative to the full 
Spitzer value. While it is difficult to constrain the globally 
averaged value of / in observed clusters, the strength of ob- 
served temperature gradients in s everal clusters suggests that 
/< 1 (IMarkevitch et alJI2003UXiang et al]i2007l) . 

There are many theoretical reasons to expect that / < 1. 
The primary uncertainty on the importance of conduction is 
the strength and order of magnetic fields in the ICM. Mag- 
netic fields act to suppress conduction across magnetic field 
lines (due to the small gyro-radius of electrons) while per- 
mitting conduction along the field lines. Observations sug- 
gest that the magnetic fields i n clusters are tan gled with 
length scales ^1-10 kpc (e.g. iTaylor et al.ir2002 h. Tangled 
magnetic fields or other magnetic phenomena tend to sup- 
press the thermal conductivity to a value roughly 10-30% 
of the fuUy Spitzer value (e . g. iMalvshkin & Kulsrudll2001l: 



ot tne tuiiy t)pitzer value (e . g. iMaivsnian & Kuisruaiizuuii: 
iNarayan & MedvedevI 120011: IChandran & MaronI |2004 . al- 
though turbulence m ay boost the effective conductivity 
jCho & Lazarianll2004l) . The amount of suppression depends 
in detail on the poorly constrained properties of the cluster 
magnetic field, and should thus be considered rather uncer- 
tain. Finally, it should be kept in mind that three-dimensional 



simulations of conduction in the presence of magn etic fields 
display new instabilities (e.g. lParrish & Stoii3l2007h that can- 
not be captured in our one-dimensional treatment where the 
conduction is by necessity isotropic. These considerations 
lead us to adopt 0.0 < f < 0.5 as a plausible range for sup- 
pression of isotropic conduction (see Table[T]). 

2.2.4. DF Heating 

A wealth of energy is stored in the orbital motions of 
galaxies within clusters. DF provides a way to transfer 
that energy to the background matter, including both the 
dark matter and the ICM, and can be a potential mecha- 
nism t hat balances the cooling flow, as first pointed out by 
lMilleil(il986) . Recen tly, DF heating has experienced a resur- 
gence in p opularity (lEl-Zant et al.l l2004bl: iKim et all 120051: 
iKimI l2007l) thanks in part to detailed cal culations of the 
efficiency of DF in a collisional medium (lOstrikerl [T999). 
which h ave been verified by controlled numerical experiments 
dSJnchez-Salcedo & Brandenburg 2001; Kim & K im 2_007J- 
These calculations showed that DF heating is stronger in col- 
lisional media (such as the ICM) when galaxies are mov- 
ing slig htly supersonically, as appears to be the case in clus- 
ters ( Falten bacher et alll2005h . where the average Mach num- 
ber of galaxies is ~ 1.3. DF will also heat the background 
dark matter halo, thereby producing a core-like inner den- 
sity profile, as opposed to the c usp implied by the NFW 
distribution (lEl-Zant et al.l l2004ah . although this appears to 
depend on the detailed properties of the accreted objects 
(BoyLan_-Kolchin & Ma 2007). 

A simple calculation suggesting that DF heating may be im- 
portant in real clusters will be presented first, followed by a 
more detailed discussion of our implementation of DF heat- 
ing. 

X-ray luminosities of clusters vary from ~ 10"*^ to ^ 10"^^ 
ergs s"', with of order 10% emitted within the rapidly cool- 
ing inner region (Fabian 1994). Assuming an average X-ray 
luminosity of lO'*'* ergs s"', this implies a loss of energy from 
X-ray radiation over a Hubbl e time of ~ 5 x 10®' ergs within 
the cooling inner region (e.g. lFabianllT99l . The energy lib- 
erated during the build-up of the massive central cD galaxy is 
a plausible source of energy that may balance these radiative 
losses. The following is a simple calculation demonstrating 
that the binding energy of the cD absorbed by the ICM is of 
the same order as the energy radiated away within the cooling 
region. 

The energy liberated during the build-up of the cD galaxy 

is 



AE= dm* =47r 



(15) 



where AE is the binding energy, (j> is the total gravitational 
potential felt by the cD, is a unit of stellar mass, and is 
the stellar density profile. We assume that cj) is static for sim- 
plicity, though of course the dark matter halo is being built-up 
at the same time as the cD. 

The cD galax y projected d ensity profile is approximated as 
a Sersic profile (iSersidll968h . in agreement with a variety of 
observations. The projected Sersic profile is not analytically 
invertible, but sufficiently accurate fitting functions exist in 
the literature. Th e de-projected stellar mass density can be 
approximated by dLima Neto et alj[r999h : 



p,(r)=As-"exp{-b„s''"), 



(16) 
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where s = r/R„ b„ w 2n-0.327, and a w 1 - 1.188/(2m) + 
0.22/(4n^). A is an integration constant chosen such that the 
integral over equals the total stellar mass M». 

For simplicity the total gravitational potential within the in- 
ner regions is assumed to be isothermal: 



0(r) = V2ln(r/reooi), 



(17) 



where Vf is the circular velocity. We are only interested in 
energy liberated within the cooling region (r < rcooi), and so 
set the potential to be zero at r = rcooi- Define 7 = r^^oox/Re- 
Combining these results yields: 

f"ln(^^)i'2-"expr-M"/"]di' 

/o i'2-"exp[-fe„i'i/«]di' 

The two integrals represent an average of the logarithmic run 
of the potential over the stellar density profile. The integrals 
are functions only of the Sersic index n, the ratio between the 
cooling radius and the effective radius 7, and the upper limit 
of integration s (taken to be the same in both integrals, which 
is reasonable but not necessary). For 1 < ^ < 3 and 4 < « < 10 
the ratio of these two integrals varies from 0.4 to 1.8, for 7=1. 
The integral scales roughly as y^. 

For a typical cD gala xy M<, ^ 10^^M(7i an d for a typical 
cluster Vc ^ 700kms"' ( Gonzalez et al. 1'2007') implying that 
/S.E ^ 10''' ergs. This energy can be transferred into both the 
background ICM and dark matter halo. If the gas were col- 
lisionless, then the fraction of energy shared between the gas 
and the dark matter would be simply proportional to their rela- 
tive mass fractions (^^ 1 /6 and ^5/6 for the gas and dark mat- 
ter respectively). However, the collisional nature of the ga s 
makes DF more efficient at transonic speeds (Os trikedl 19991) : 
at Mach numbers near unity the efficiency of DF heating on 
the gas is roughly twice that of the collisionless dark matter. 
This implies that ~ 30% of the energy released in transferred 
to the ICM and the rest to the dark matter. The cD binding 
energy thus absorbed by the ICM is comparable to the energy 
radiated away within the cooling inner region over a Hubble 
time. 

There were many simplifications made in this calculation. 
One complication is that observations indicate that a large 
fraction of massive ( M^ > 10"^M(y,) gal axies were already 
in place by z ~ 1 (e.g. 'B undv e"tlLll2005l: iBorch et al.ll2006i 
iFontana et al.l l2006: Brow n et al.ll2007l) ."and hence the energy 
liberated from the build-up of cDs may not provide much en- 
ergy to the ICM at z < 1. However, the outer envelope of 
the cD may be growing substanti ally at z < 1 from the shred- 
ding of satellite galaxies (see e.g.[Ostriker & Hausmanll 19771: 
IConrov et al.ll2007l:lPurcell et al.ll2007l for a discussion), sug- 
gesting that energy released from the build-up of the cD is 
available at late times. Moreover, the strong small-scale clus- 
tering of luminous red galaxies indicates that there are many 
massive galaxies orbiting near the cluster center, plausibly 
transferring their orbital energy to the ICM and background 
halo via DF (Masjedi et al."2006). 

These calculations suggest that DF heating may be an im- 
portant heating mechanism. While the simulations below 
address DF heating in more detail, detailed high-resolution 
three-dimensional simulations are required to fully address 
the importance of DF heating in the ICM of clusters (see e.g. 
|Kim 2007). Moreover, detailed resolution studies are required 
to address whether or not current cosmological hydrodynamic 
simulations are adequately resolving DF (see discussion in 
iFaltenbacher et al.ll2005HNaab et al.ll2007h . 



Our implementat ion of DF heating closely follows that 
of iKim et al.l (l2005l) : the reader is referred to that work for 
more details. Note that we do not track orbits explicitly nor 
do we attempt to resolve the actual DF wake. We take a 
more approximate, parameterized approach to the spherically- 
averaged heating rate due to DF. 

The heat flux due to DF can be described as: 



rDF = ngal(-FDF-v) 



A^p.G^MUlIM) 



ngai(r), (19) 



where Ai = v/c, is the Mach number, n^.^\ is the number den- 
sity of galaxies, M^^j is the average squared total mass of the 
galaxies, and c, is the isothermal sound speed. In what follows 
we will argue for reasonable values for each of the parameters 
in Equation[T9l and will then incorporate the uncertainties and 
possible ranges in all of these parameters into the single free 
parameter d. The angular brackets indicate an average over 
the velocity distribution function. 

Note that Mgai includes both the stellar and dark matter of 
the satellite galaxies, and that since DF is proportional to the 
average of the galaxy mass squared, more massive galaxies 

are given greater weight. The quantity M^^j is estimated as 
follows. We assume that the dark matter-to-stellar mass ra- 
tio is constant for satellites, and estimate it by taking the total 
stella r mass of satellites in c lusters to be 1 % of the total cluster 
mass jGonzalez et alj20"07D and the t otal amount of da rk mat- 
ter mass in substructures to be 10% dGao et al.ll2004l) . where 
the substructures are assumed to be the likely locations of 
satellite galaxies (Conroy et al. 2006J. These numbers imply 
an average dark matter-to-stellar mass ratio of 10 for satellite 

galaxies. If the total mass is Mgai, then M^^, = MI^.^^ + M\^ = 
lOlMstar- The stellar mass function provides an estimate of 
Mjtaj and thus of M^^j. We adopt the global mass function 
of lBeUetalJ (l2003h where a = 1.1 and log(M*) = 10.9Mo 
are the best-fit Schechter parameters. The Schechter parame- 
ters for the luminosity function within clusters does not differ 
strongly from the global value dHansen et al.l 1200 7). and so 
we adopt the global values for this calculation. Integrating 
the stellar mass function to Mjtar = 10^ Mq leads us to adopt 
M^jjj = (10" Mq)^ as the fiducial value. The uncertainty on 
this quantity is explored via the tunable parameter d discussed 
above. 
By definition we have: 



{I/M): 



J I/Mf(v)dv 
J mdy ' 



(20) 



where / is assumed to be Maxwellian (IFaltenbacher et all 
l2005h : 



/(v) = 



(2^a2)3/2 



(21) 



where cr, is the radial velocity dispersion of galaxies (assumed 
for simplicity to be independent of radius), and A^gai is the total 
number of galaxies within the cluster. In the following we 
take cr,. = 1000 km s"' and A'gai = 500. These parameters are 
fixed throughout; any variation of these physical quantities 
is incorporated into the d parameter. Note that DF heating 
depends on the state variables only through 1/ M., and pg, 
and that Vof « pg at fixed temperature. / is the efficiency 
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TABLE 1 
Summary of Free Parameters 







Plausible 


Parameter 


Comment 


Range 


Z/Zo 


Metallicity of the gas (Eqn. I12t 


0.1-0.6 


/ 


Fraction of Spitzer conductivity (Eqn. I14t 


0.0-0.5 


d 


Normalization of DF heating (Eqn. |19[ 


0.1-10.0 


a 


Initial fraction of relativistic pressure (Eqn. I23t 


0.0-0.3 



iln(l-A^-2) + ln(vt/r„i„) > 1 



l-M 



-M 



M < 1 



As is common practice, we set vt = rmax, where r„iax and r^in 
are in this problem taken to denote the size of the cluster and 
the satellite galaxies, respectively. The factor ln(vr / rmin) thus 
plausibly varies from ~ 4 - 10. In what follows we set this 
factor to 6 and note that our results are insensitive to this par- 
ticular value. Note that DF heating is more efficient at sub- 
sonic speeds than one might naively expect because the rele- 
vant efficiency is {I/A4), which, though maximal at = 1, 
decreases by only 50% at Al « 0.5 and decreases only weakly 
at > 1 . This implies that the feedback provided by the 
Mach number (in the sense that colder/hotter systems will 
have more/less efficient DF heating) is weaker than one might 
have expected. 

The number density of galaxies, «gai(?")> is taken to be 
a modified King p rofile with the parameters adopted from 
iGirardi et aUdlOOSh : 

«gai(r) = «gai(0)[l + (r/r,)2]-'-2, (22) 

where the core radius = 50 kpc. The central density is set by 
requiring that «gai(?") integrate to the total number of galaxies, 
A^gai- Note that n„i^\{r) does not include the central cD galaxy. 
While most cD galaxies are near the center of the halo as de- 
fined by X-ray imaging ( Lin & Mohrll2004b . the cD may be 
oscillating abo ut the gravitational center with an amplitude 
of several kpc (Ivan den Bosch et ani2005h . Such oscillations 
may provide additional heating in the central regions, but we 
do not include them herein. 



2.3. Relativistic Pressure 

Relativistic pressure in the form of cosmic rays could be 
an important dynamical component of the ICM. It has long 
been known that astrophysical shocks e fficiently accelerate 
cosmic rays (B landford & Ostrikeiill978 ). Numerical simu- 
lations suggest that cosmic rays accelerated in shocks could 
account for ^ 10% of t he total cluster pressure (iMiniati et al.l 
l200lt iRvu et alJ l2003h . Cosmic rays may also be gener- 
ated in buoyantly ris ing bubbles generated by AGN activ- 
ity (e.g. 'E nBUnll2003h . Observations suggest that a relativis- 
tic component may constitute seve ral tens of pe rcent of Jhe 
total energy density of clusters (jP frommer & EnBlin 2004; 



Sanders etal]l2005 



20071 IWerner et al 



iDunn & Fabian 2006; Sa nders & Fabian 
|2007; Nakaret al. 2007). In addition, it 
has been shown that a relativi stic pressu re component can 
suppress thermal instabilities (ICenI 120051 see also the ap- 
pendix herein). 

There have been several recent attempts to study the ef- 
fects of cosmic-rays in the ICM numerically dPfrommer et alJ 



2007blal: iPfrommedllOOTl: TGuo & OhilIOOTh . The model of 
Guo & OhI (I2007h . which includes both conduction and cos- 



mic rays injected into the ICM from AGN-induced bubbles, 
is able to reproduce the temperature and density profiles of 
observed clusters. Perhaps the most attractive feature of their 
model is that their results do not require fine tuning of the var- 
ious adjustable parameters, including the amount of thermal 
conductivity. It would be interesting to know whether and to 
what extent their results rely on the relativistic pressure pro- 
vided by the CRs, or whether it is due primarily to the ener- 
getics associated with the bubbles. In the present work we test 
the former hypothesis. 

In the following we take a simplified approach when ex- 
ploring the effects of a relativistic component. The relativistic 
pressure (P, ) is assumed to be a fixed fraction of the total pres- 
sure initially: 

Pr = aP,o„ (23) 

and we further assume that the relativistic component is per- 
fectly dynamically coupled to the gas (i.e. there is no mo- 
tion or diffusion of one component relative to another), and 
evolves adiabatically (i.e. Pr cx pj' ). These requirements lead 
to a fourth fluid equation: 



dp 

_Z=_v.(vP,) + (l- 



(24) 



which is simply a statement of energy conservation. In the 
above equation 7,. = 4/3 and we have made use of the follow- 
ing equation of state: 

^ :Pr. (25) 



1 



where is the energy density associated with the relativistic 
fluid. Note that relativistic pressure also enters into the hydro- 
dynamical equations (Equations [T]—|3]l by contributing to the 
total pressure. 

These simple assumptions are motivated by CR creation 
and evolution in real clusters. It may be the case that CRs are 
created predominately in the cluster center via AGN activity, 
or it may be that the accretion shock at the cluster outskirts 
are the predominate source of CRs. CRs are capable of dif- 
fusing out of the cluster, but CR loss via diffusion may be ap- 
proximately balanced by the creation of new CRs throughout 
the cluster. Our simplified treatment is meant to demonstrate 
the potential importance of CRs generally (insofar as they 
provide relativistic pressure support); it wil l be the task of 
more sophisticated model i ng efforts (see e.g . iPfrommer et alj 
'2007b a; Pfrommer "2^ |Guo & Oh 2007) and, ultimately, 
observations to refine our knowledge of CR production and 
evolution. Note that the energetics associated with CR pro- 
duction (whether in e.g. shocks or AGN-related bubbles) is 
not explored herein. 

2.4. Initial Equilibria 

Every simulated cluster is set up initially in hydrostatic 
equilibrium. However, in general the clusters are not set up 
in thermal balance. 

Observations indicate that there are two classes of intra- 
cluster media, t he cooling-core (CC) and non-cooling-core 
(NCC) clusters (ISanderson et al.ll2006l) . Clusters in the for- 
mer class show a rapidly declining temperature profile with 
decreasing radius at radii less than ~ 100 kpc, while the lat- 
ter are approximately isothermal within the same physical re- 
gion. Furthermore, CC clusters have factors of ~ 2 - 3 metal- 
licity gradients within ^100 kpc while NCC clusters do not 
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dPe Grand! & MolendilllOOl . For our analysis we assume 
that neither NCC nor CC clusters have an appreciable metal- 
licity gradient, and that the metallicity is (1 /3)Z0. We have 
selected several cluster runs at random and re-simulated them 
with metallicity varying from primordial to solar composition 
and find no qualitative change in our results. 

Both of these classes are discussed in the sections that fol- 
low. NCC clusters as approximated as initially isothermal 
while the CC clusters are assumed to have an initial tem pera- 
ture profile that is described by ( IVikhhnin et al.ll2006bh : 

1.15(jc/0.045)2+0.5 1 



T{r)/Tcc ■■ 



(26) 



(x/0.045)2+l l + (x/0.75)2' 

where x = r/ rsoo. This temperature profile peaks at Tqc where 
r ^ O.lrsoo and drops by a factor of ~ 2-3 toward the center. 
Figure |2](/7oftom panel) provides an example of a CC cluster 
temperature profile. 

In each run the cumulative gas fraction at 500 kpc is 
fixed at /e = 0.1, in agreement with recen t observations 
(ISanderson et al.ll2003HVikhlinin et alJl2006bl) . The gas frac- 
tion is less than the univers al baryon fraction of ft = 0.16 
for several reasons (see e.g. iKravtsov et alll2005L for a dis- 
cussion). The conversion of gas into st ars results in a lower 
gas fraction by several tens of percent (iFukugita et al.|[T998l : 
iLin et alj|2003h . Both adiabatic effects and heating processes 
furthe r decrease the gas fraction within ^ 500 kpc (Bode et al. 
120071) . Fixing the gas fraction implies that simulated clusters 
with larger central electron de nsities are embedded in larger 
dark matter halos, as observed dVikhlinin et alj|2006bh . 

With the above simplifications, each simulation is fully 
specified by the initial central electron density and the tem- 
perature profile. These two parameters, along with the as- 
sumption of hydrostatic equilibrium, then determines the full 
density, temperature, and pressure profiles, and the adopted 
gas mass fraction then determines the mass of the static dark 
matter halo. The contribution of the cD galaxy to the gravi- 
tational potential is fixed initially and increases in proportion 
to the amount of mass that flows through the inner bound- 
ary condition (see !j2.1l ). In practice we iteratively solve for 
Mo such that hydrostatic balance is achieved. For the central 
densities explored below, the NFW normalization Mq ranges 
from ^ (2-7) X 10 '"^Mq. AUowing Mq to vary by fixing fg as 
opposed to allowing fg to vary by fixing Mq has a negligible 
effect on our conclusions. 

Observationally the temperature of clusters varies from ~ 
1-10 ke V, and is strongly cor related with the total cluster 
mass (e.g. Vikh linin et al. !'2006b). The data also cover a wide 
range in central electron densities, roughly 10"^^ < «(,,o < 
10""^ cm"^. We pay more attention to how cooling and heat- 
ing depends on density because the cooling and heating func- 
tions are more sensitive to density than temperature. Below 
we include "observed" central electron densities for compari- 
son with our initial conditions, which a re estimated from two 
sources. IZakamska & NaravanI (l2003h provide central den- 
sities from hydrostatic fits to obse rved temperature and den- 
sity profiles of 10 clusters. Vikhli nin et alj ( l2006bh provide 
detailed parametric fits to the electron density profiles of 13 
clusters (two clusters are in both samples and are only counted 
once herein). For the latter sample we quote as the central 
density the density at 10 kpc because in many cases the den- 
sity profile at smaller scales is enhanced due to cooling. Since 
we are interested in "initial" central electron densities, we es- 
sentially mask out this inner cooling region when comparing 
to the initial electron densities in our simulations. In the simu- 



t = 0.0 Gyr 




t = 0.0 Gyr 




1.5 2.0 2.5 3.0 
logio(r [kpc]) 

Fig. 2. — Evolution of the temperature profile of clusters with radiative 
cooling and no heating for n^ Q = 0.02 cm"^^. As expected, a runaway cooling 
flow develops and the core collapses within a cooling time. Top Panel: An 
initially NCC cluster. Bottom Panel: An initially CC cluster. The lines are 
spaced equally in time. 



lations the initial central electron densities are approximately 
constant within the inner 10 kpc; we thus consider this sim- 
plification appropriate. 

2.5. Numerical Setup 

We utilize the Lax-Wendroff method dPress et al.ll 19921) to 
integrate the fluid equations in their one-dimensional form. 
The spatial grid is logarithmic and has = 400 elements with 
range 1 < r < 500 kpc. We have doubled the number of grid 
elements for several simulations and extended the spatial grid 
to 1000 kpc and find no change in the results. The time-step 
is determined by the Courant condition: 



AfcFL = 0.5 m/n 



\\i + Cs.i\ 



(27) 



where ; = [0,A^] and , is the local isothermal sound speed. 

Conduction is i mplemented with a fully implicit algorithm 
(iPress et al.lll992l) . The pressure, density, and temperature in 
the ghost cells are linearly extrapolated from the active zones, 
both at the inner and outer boundaries. At the inner boundary, 
the velocity is also a linear extrapolation from the active zones 
unless the velocity in the ghost cell is positive (indicating out- 
flow) in which case the velocity is set to zero. The velocity in 
the ghost cell at the outer boundary is always zero. 

3. NUMERICAL EXPERIMENTS 

This section presents the results of a series of numerical ex- 
periments. The various physical processes discussed in ^ 12. 21 
are discussed in a variety of combinations in order to elucidate 
their various effects. For each process we discuss the sensi- 
tivity of the final state to the free parameters summarized in 
Table [U 
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4 6 
time [GyrJ 

Fig. 3. — Evolution of the central temperature for simulations that include 
both radiation and DF heating. Each panel shows the evolution of clusters 
with a variety of initial central densities (labeled in the figure in units of 
cm"'). The normalization of DF heating here is fixed at d = 1.0 (cf. Equation 
\l9\ . Top Panel: Initially NCC clusters. Bottom Panel: Initially CC clus- 
ters. The evolution of CC clusters is qualitatively similar to NCC clusters in 
the sense that thermal balance can only be maintained for a narrow range of 
central electron densities. 



3.1. Radiation Only 

As discussed in the Introduction, in the absence of any heat- 
ing mechanisms the ICM will cool catastrophically. The time 
required to radiate away all of the internal energy of the gas is 
defined as the cooling time': 



^cool 



~ A 



:5.0 



0.01 cm" 



—) 

6keV J 



0.5 



Gyr. 



(28) 



As can be seen from Equation|28] the central cooling time for 
typical clusters is less than a Hubble time. 

Figure |2] shows the evolution of the ICM for two clusters 
with radiative cooling and no heating. It is clear that the 
clusters cool catastrophically in the absence of any heating. 
The top panel shows the evolution of an initially isothermal 
cluster with T = 6 keV and an initial central electron den- 
sity «(..() = 0.02 cm"^; the bottom panel shows the evolution of 
an initially CC cluster with the same central electron density. 
The cooling times for the NCC and CC clusters shown in the 
figure are 2.4 and 1.7 Gyr, respectively, based on Equation 
l28l and adequately captures the actual time for catastrophic 
collapse. 

While there are clusters with temperature profiles simi- 
lar to that shown in the top panel of Figure |2] it is impor- 
tant to note that not only are these temperature profiles tran- 

' Note that this cooling time differs from the conventional cooling time 
by a factor of because the conventional cooling time is defined as the 

time required to cool isobarically whereas the cooling time in Equation |28| is 
simply the time required to radiate away all of the thermal energy. 



sient (in the sense that the cluster continues to cool on rapid 
timescales), but the density profile of this radiation-only clus- 
ter (not shown) is nowhere observed in nature. Hence the 
X-ray luminosities, which scale as the gas density squared, 
would be out of the range of those observed in nature. The 
conclusion that observed clusters must be periodically heated 
is in accord with a large body of previous work. In partic- 
ular, X-ray spectroscopy has firmly established that the ICM 
is not cooling much below^ 1 /3 of the ambien t temperature 
dPeterson et alJl200Ul200llTamura et alll2001h . 

The following sections discuss possible heating mecha- 
nisms that can forestall the cooling catastrophe visible in Fig- 
ured 

3.2. Radiation & DF Heating 

As discussed in !j2.2.4l there is an enormous source of en- 
ergy in the orbital motions of satellite galaxies. While satellite 
galaxies can transfer their orbital energy to the ICM through a 
variety of mechanisms, we will focus only on DF heating for 
simplicity. Note however, that other mechanisms capable of 
transferring the orbital energy of satellites to the background 
gas are expected to scale in a similar way as DF heating (i.e. 
as pg at fixed temperature). The discussion that follows thus 
roughly encompasses a variety of heating mechanisms related 
to the motions of satellite galaxies. 

Figure [3] shows the evolution of the ICM central tempera- 
ture, Tc, as a function of time for initially NCC (i.e. isother- 
mal; top panel) and CC clusters (bottom panel). Each panel 
shows the central temperature evolution for a range of initial 
central gas densities. The normalization of DF heating shown 
in the figure is set to li = 1 .0. 

Several trends are apparent. First, it is clear that one can 
find a particular set of parameters that leads to approximate 
thermal balance (neither a heating nor cooling catastrophe) 
over a Hubble time. While not shown, we note that in the 
runs where the central temperature does not change appre- 
ciably, neither do the density nor temperature profiles. This 
was first demonstrated by Kim et al. (.2005 ). In fac t, our fidu- 
cial va lues for DF heating are the same as those in lKim et alj 
(l2005h . The equilibrium NCC model has Tij^cc = 6 keV and 
«e,o = 6 X lO""' cm~ ^ which is exactl y the set of equilibrium 
parameters found in lKim et alj ( l2005h . thereby confirming the 
results of that work. 

It is also apparent that DF heating as the sole heating mech- 
anism leads to serious thermal imbalance for all but a very 
narr ow ran ge of electron densities, as was also demonstrated 
bv lKim et alj ( l2005h . This imbalance arises because at fixed 
temperature the cooling scales as while DF heating scales 
as Pg. Thus at fixed temperature there can be only one initial 
density that is in thermal balance, and that solution will be 
unstable. 

In order to more fully elucidate the thermal properties of 
the ICM, Figure |4] shows the long-term thermal evolution of 
the ICM as a function of initial central density, n^.o and nor- 
malization of DF heating, d, for NCC (left panel) and CC 
(right panel) clusters for a fixed initial temperature profile. In 
this figure the slope of the central temperature as a function 
of time is color-coded according to the sign and steepness of 
the slope. From Figure [3] it is clear that Tc for most clusters 
evolves at least quasi-linearly, and hence the slope is a reason- 
able metric for the evolution of the ICM. 

Figured demonstrates that there is only a narrow range of 
thermal balance in the parameter space of d and n^.o. This 
highlights the inability of DF heating operating alone to bal- 
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Fig. 4. — Long-term thermal evolution of inti' aclus ter media subjected to both DF heating and radiative cooling, as a function of initial central density, n^ o 
and normalization of DF heating, d (cf. Equation ll9t . Both NCC clusters with Tj ncc = 6 keV (left panel) and CC clusters with Tj cc = 6 keV (right panel) are 
shown. The color indicates the slope of central temperature vs. time. Red indicates the region of parameter space where the cluster gas catastrophically heats 
and blue where the gas catastrophically cools (i.e. where the slope is positive or negative). The solid blaclc lines indicate the region of parameter space where the 
slope is equal to ±0. 1 . Note that while there is a domain where heating stably balances cooling, sensitive fine- tuni ng is required for a real cluster to remain in 
this domain. For comparison, the distribution of observed central electron densities is included at the top (see ^2.4) and the estimated location of several Abell 
clusters are included in the d—n^fl plane (circles). 
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time [GyrJ 



Fig. 5. — Evolution of the central temperature for simulations that include 
both radiation and conduction. The conduction normalization is / = 0.1 (cf. 
Equation lMt . Top Panel: Initially NCC clusters with Tj j^cc - 6 keV. Bottom 
Panel: CC clusters with 7} cc = 6 keV. 



ance radiative losses in the ICM. The normalization of DF 
heating that is required to balance radiative losses increases 
in proportion to the central electron density id cx n^.o), which 
makes the DF heating term (Equation [T9]l effectively scale as 



p^. This is not surprising because the cooling function scales 
as and so thermal balance can be achieved for a range of 
densities if the effective DF heating scales in the same way as 
cooling. 

Nature may in fact provide this required scaling {d oc tie o). 
For example, the number of galaxie s per cluster is correlated 
with the cluster mass as A^gai cx M*' ** (Lin et al.'2004'), with ap- 
parently little change at higher redshift (Lin et al. 2006), and 
the central electron density scales rou ghly as the cluster mass: 
We (X M (though with larg e scatter: IVikhlinin et al.l l2006bl : 
IZakamska & Narayarill2003h . Since we have incorporated the 
variation of all free parameters in Equation [19] into the single 
parameter d, it follows that d cx A^gai- These additional scal- 

ings may thus in fact produce d oc q with [3^1. 

In order to investigate this explicitly, we have plotted the 
approximate locations of several Abell clusters in the cZ-n^o 
plane in Figure 21 w here the centr al electron density is esti- 
mated from Vikhlinin et alJ (l2006bl) and the number of galax- 
ies, Nm\, is provided by the luminosity functions of iLin et aH 



(|2004|) . We have only included the dependence of d on A^g^i 
when estimating the different DF efficiencies for these clus- 
ters since this scaling is the dominant one. In other words, 
for these clusters we adopt our fiducial values relevant for 
DF heating except for A^g^i, which we take from the literature. 
Moreover, the intracluster media of these clusters span a range 
in temperature (from ^2-9 keV; Lin & Mohr 2004). Since 
the region of thermal balance in Figure|4]is only a weak func- 
tion of temperature, we include all the clusters on the same 
plot for simplicity. A trend of increasing d with increasing n^, o 
is evident, but the normalization is a factor of ^ 3 lower than 
the region where DF heating balances cooling. While a factor 
of three is probably within the uncertainty of our estimates of 
d for these clusters, what is more important is the narrowness 
of the region. Our estimates of d would not only have to be 
systematically off by a factor of ^ 3, but there would have to 
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Fig. 6. — Same as Figure |4] except now only thermal conduction and ra- 
diative cooling are considered, and the clusters are initially all NCC. The 
normalization of the conductivity, / (cf. Equation ll4t . is varied as a function 
of the initial central electron density tif Q. The initial gas temperatures are 5 
keV (top panel) and 8 keV (bottom panel). 



be rather extreme fine-tuning, especially at high n^. o, for DF 
heating to be generically be able to balance radiative losses. 
Moreover, even if clusters initially fell in this narrow range, 
DF heating would still respond to perturbations in q only 
linearly, while cooling would respond quadratically, and thus 
it is unlikely that clusters would remain in that narrow region 
for long. 

The conclusion drawn fro m this section is st raightforward 
and echoes the conclusion in ' Kim et alj (l2005h : DF heating, 
though an important reservoir of energy, cannot be the sole 
heating mechanism operating to offset radiative losses in the 
ICM. This is not because the energy available is insufficient 
but rather because there is only a narrow range of parameter 
space where DF heating can stably balance radiative cooling. 

3.3. Radiation & Conduction 

Thermal conduction can in principle transport heat from the 
abundant thermal reservoir of the outer cluster gas to the inner 
cooling parts. Invoking conduction as a means to balance ra- 
diative losses in the ICM appears howe ver to require a fair de- 
gree of fine tuning of the conductivity dZakamska & NaravanI 



iKim & NaravanI l2003i iPope et all 120061: IGuoXoE 
As demonstrated below, this problem is particularly 
acute in CC clusters where the observed temperature profile 
drops precipitously in the inner parts — too little conduction 
leads to a cooling catastrophe, while too much conduction 
tends to produce an isothermal core. 

These issues are demonstrated graphically in Figure |5] 
where the evolution of the central temperature is plotted for 
both NCC {top panel) and CC {bottom panel) clusters, for a 
conductivity normalization / = 0.1. Notice that the evolution 
is qualitatively different for NCC and CC clusters. 

For NCC clusters, which are initially isothermal, only clus- 
ters with very low gas densities are marginally stable, while 
progressively more dense clusters quickly runaway. These 
trends can be understood with the aid of the conduction time: 

2 /^, ,,\ 2.5 



04 

T 



0.01 cm-3 J VlOOkpc 



ekevy 



Gyr, 



(29) 

where / is the conductivity normalization and A is the length 
scale over which temperature gradients are washed out within 
a conduction time (the inner region where the cooling time 
is shorter th an a Hubble time is generally of order 100 kpc; 
ISandersolTe t al. 2006). It is clear that larger densities lead to 
longer conduction times^, and that for / = 0. 1, densities much 
greater than n^ o ^ 10"^ cm"-' result in conduction times ap- 
proaching the Hubble time, and much longer than the cooling 
time. For these high densities, conduction cannot forestall 
runaway cooling, as seen in Figure |5] Thermal conduction 
alone cannot therefore stably offset radiative losses in the in- 
tracluster media of NCC clusters. 

The situation is somewhat more complex in CC clusters be- 
cause there the initial temperature profile drops by a factor 
of 2-3 within the inner ^100 kpc. The bottom panel 
of Figure |5] shows the central temperature evolution for these 
clusters. For clusters with short conduction times, the inner 
(cooler) region isothermalizes to the temperature of the outer 
cluster region 6 keV for these runs) before radiative effects 
become important. Once isothermal, the clusters then evolve 
more like the NCC clusters {top panel. Figure |5]l. CC clus- 
ters that are too dense have conduction times that are longer 
than the cooling times, and they thus cool catastrophically. 
These trends persist for the full range of conductivities ex- 
plored herein (0 < / < 1). 

The trends evident for NCC clusters in the top panel of Fig- 
ure |5] are shown for a wide range in conductivity normaliza- 
tions and central electron densities in Figure |6] We have sim- 
ulated NCC clusters with an initial temperature of 5 keV {top 
panel) and 8 keV {bottom panel) in order to demonstrate ex- 
plicitly that our results are relatively insensitive to the ICM 
temperature. Unlike Figure |4] NCC clusters can only cool 
or remain in thermal balance, since conduction acting on an 
isothermal gas cannot produce a steadily increasing central 
temperature with time. It is clear that increasing the conduc- 
tivity results in stable intracluster media over a wider range 
of initial central electron densities, but the range where cool- 
ing balances heating increases only mildly as the conduction 
normalization increases from / ^ 0.2 to / ~ 1 .0. 

The DF heating and conduction fine tuning problems are 
thus somewhat different. In the former case, too much DF 
heating will result in runaway heating, while in the latter, too 

^ The conduction time depends on density because, while conduction 
transports energy at fixed density, energy scales as n^T, so that the conduction 
time scales as n,. 
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Fig. 7. — Evolution of the central temperature for simulations that include 
radiation, DF heating, and conduction. The DF heating and conductivity nor- 
malizations are d = i.O and / = 0. 1 , respectively. Top Panel: Initially NCC 
clusters. Bottom Panel: Initially CC clusters. 



much conduction simply isothermalizes the core. This is not 
acceptable because the majority of clusters s how a factor of ^ 
2-3 temperature drop within their cores (e.g. lSanderson et aTl 
I2006t IVildilinin et alJl2006bh . The implications are however 
the same: neither mechanism is generically able to reproduce 
observed temperature profiles. 

It thus appears that previously proposed models for the 
ICM that construct clusters initially in thermal b alance with 
conduction as the only heating source jZ akamska & NaravanI 
I200I iKim & NaravanI l2003atlGuo"&Ohi i2007) constitute a 
very narrow regime of parameter space and should thus not 
be considered realistic solutions to the cooling flow problem, 
unless a physical explana tion for the fine-t uning is provided 
(see also the discussion in lGuo & 0^120071) . 

3.4. Radiation, DF Heating, & Conduction 

In this section the thermal balance of the ICM is explored 
when both DF heating and conduction are operating as heat- 
ing mechanisms. 

Figure I2] plots the evolution of the central temperature for 
NCC and CC clusters with normalization of DF and conduc- 
tion set to d = 1.0 and / = 0.1. Comparing this figure to 
Figures |5] and [3] highlights the stabilizing effects of the com- 
bination of DF heating and conduction. Note however that 
the evolved CC clusters are plagued by the same issues dis- 
cussed in ! j3.3l namely that while the combination of DF and 
conduction can maintain CC clusters in thermal balance, the 
evolved clusters do not preserve the factor of ~ 2 - 3 drop in 
the temperature profile characteristic of observed CC clusters. 
In other words, these CC clusters that do not catastrophically 
cool instead turn into NCC clusters, and thus these runs can- 
not explain the existence of observed CC clusters. 
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Fig. 8. — Same as Figure |4] now with DF heating, thermal conduction, 
and radiative cooling, for NCC clusters with Tj ^cc = 6 keV. Top Panel: NCC 
clusters with the conduction normalization set to / = 0.1. Bottom Panel: 
NCC clusters with / = 0.5. Note that even a small amount of conduction 
makes the clusters much more thermally stable than DF heating alone over a 
much wider range of central gas densities. As in Figure |4]the approximate 
locations of several Abell clusters are included (circles). 



The increased thermal balance when combining DF heat- 
ing with conduction is manifest in Figure [8] which shows the 
evolution of the central temperature of the ICM for initially 
NCC clusters with Ti ffcc = 6 keV in the parameter space of 
initial central electron density and normalization of DF heat- 
ing. This figure shows the effect of including conduction with 
normalization / = 0. 1 {top panel) and / = 0.5 (bottom panel). 

Upon comparing this figure to Figures |4] and |6] it is clear 
that the combination of DF heating and conduction makes the 
ICM of NCC clusters considerably more thermally balanced 
than either mechanism acting alone. Indeed, for initial central 
electron densities n^ o S 02 cm"^ the ICM is thermally bal- 
anced for any plausible value of the DF heating normalization, 
and for reasonable, not fine-tuned, choices of conductivity. 

However, it is still the case that for ne o > 0.02 cm""*, no 
plausible amount of DF heating can offset the cooling catas- 
trophe. Increasing the amount of conduction to / = 1.0 ex- 
pands the zone of thermal balance to the right in Figure [8] 
only marginally. While some observed NCC clusters have a 
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central density lying within the region of thermal balance in 
Figure |8] it is clear that for many observed clusters processes 
in addition to DF heating and conduction are at work in order 
to effectively balance radiative cooling. 

These conclusions are rather insensitive to the initial tem- 
perature of the NCC and CC clusters. This might be expected 
because the cooling function scales as r'/^ while DF heating 
scales at r"'/2+'' with p ^0 dKim et al.ir2005h . For exam- 
ple, decreasing the initial temperature in Figure |8]by a factor 
of two results in somewhat more thermally balanced clusters 
in the high density regime. However, the increased zone of 
thermal balance is narrow at high density and hence rather 
fine tuning would be required for an observed cluster to re- 
main there. Moreover, observed clusters with larger central 
densities have higher rather than lower temperatures, and so 
the problem of generating equilibrium clusters is only exacer- 
bated if the observed temperature-density relation is consid- 
ered. 

In sum, while it appears that the combination of DF and 
conduction can produce intracluster media in thermal balance 
for a modest range of central electron densities for NCC clus- 
ters, neither NCC clusters with moderately high central elec- 
tron densities nor CC clusters of any density can be kept in 
thermal balance for cosmological timescales with these mech- 
anisms alone. 

3.5. The Impact of a Relativistic Fluid 

A relativistic fluid may be present in observed clusters, per- 
haps in the form of cosmic rays ([ Pfrommer & EnBlin 2004; 
Sanders et al. 2005; Dunn & Fabian 2006; Sanders & Fabian 
2007; Werner et al. 2007; Nakar et al. 2007). We show in the 
appendix that the addition of a small amount of relativistic 
pressure increases the thermal stability of a gas already in 
thermal balance (i.e. the growth of thermal perturbations is 
damped for intracluster media where the net heating of a par- 
cel of gas is zero). In this section we explore the impact of a 
relativistic fluid on the t hermal balance of the ICM. 

As discussed in ^12.31 a relativistic fluid is included in our 
simulations by requiring that its pressure be a fixed fraction of 
the thermal pressure initially. The relativistic fluid is assumed 
to be massless and perfectly dynamically coupled to the gas. 
Its subsequent adiabatic evolution is then governed by energy 
conservation. Recall that in our treatment the relativistic fluid 
provides no additional thermal energy to the system. The pur- 
pose of this experiment is thus to ascertain the potential im- 
portance of a relativistic fluid on the hydrodynamics of the 
ICM, not the thermodynamics. 

We have run a series of simulations that span the range of 
parameter space shown in Figure |4] i.e. for a range in ini- 
tial central electron densities and DF heating normalizations. 
In addition to varying these two variables, we have also var- 
ied the initial fraction of relativistic pressure, a, from 10% to 
40% of the total pressure, which is a fraction con sistent with 
observations ([P frommer & EnBlin 2004; Sanders et al.ll2005l; 
Dunn & Fabian 2006; Sa nders & Fabianii2007c i Werner et al.l 
2007; Nakar et al. 2007). We find that flie addition of rela- 
tivistic pressure has a completely negligible effect on the ther- 
mal balance of the gas for the entire range of parameter space 
we have explored. 

At first glance this is somewhat surprising since a thermal 
stability analysis indicates that a small amount of relativistic 
pressure should increase markedly the stability of a gas al- 
ready in thermal balance (Cen 2005, see also f|A]i. The key 
difference is that stability analyses assume the gas to be in 



thermal balance, while the simulations we have run are al- 
most never set initial l y in t hermal balance. In the former 
case the classic iFieldl (Il965h instability (described in the In- 
troduction) is suppressed because slightly over-dense fluid el- 
ements, which thus cool rapidly and hence lose thermal pres- 
sure support, contract less than they would otherwise thanks 
to the additional pressure support provided by the relativistic 
fluid. However, in the latter case, where the gas in not in ther- 
mal balance, the additional pressure support provided by the 
relativistic fluid is immaterial because the thermal runaway 
is caused by a more serious imbalance between heating and 
cooling, which the relativistic fluid, providing only pressure 
support and not energy, cannot remedy. 

4. DISCUSSION 

The results of the preceding section imply that both thermal 
conduction and DF heating are potentially significant heating 
processes in the ICM. In fact, for large regions of parame- 
ter space, these processes provide more than enough energy 
to offset radiative losses. However, we have found that they 
can neither produce nor maintain intracluster media in thermal 
balance over cosmological time-scales. In particular, these 
processes tend to generate either runaway cooling or run- 
away heating, in both cases producing temperature and den- 
sity profiles of the ICM nowhere near those observed in nature 
([Sanderson et al. 2006; Vikhlinin et al. 2006b). These con- 
clusions expand upon those of iBrighenti & Mathewsl (1200^ 
who found that none of a large range of possible steady-state 
heating mechanisms could generically reproduce the observed 
properties of the ICM. For the case of DF heating, these re- 
sults are driven largely by the fact that radiation losses scale as 
the gas density squared, while DF heating scales linearly with 
the density. Since observed intracluster media exhibit a wide 
range of densities, fine-tuning is required to prevent thermal 
runaways when DF is the dominant heating mechanism. 

These conclusions have several implications, two of which 
we now discuss in detail. First, in light of our results from 
one dimensional simulations, we discuss the continued neces- 
sity of high-resolution, fully three-dimensional, cosmologi- 
cally embedded hydrodynamic simulations for understanding 
the thermodynamics of the ICM. Second, we discuss a heat- 
ing mechanism that may plausibly generate and maintain the 
properties of observed intracluster media, since none of the 
mechanisms explored herein appear capable of doing so. 

4.1. The Cosmological Context 

Our approach has been to simulate in one dimension the 
evolution of the ICM when subjected to various heating mech- 
anisms and radiative cooling in static, isolated clusters, where 
only the ICM was evolved. In the real Universe, however, 
clusters do not appear to evolve in isolation. Indeed, cosmo- 
logical cold dark matter A^-body simulations demonstrate that 
cluster-sized dark matter halos (M > 10'"* M©) continually ac- 
crete additional halos with a range of masses, many of which 
likely host galaxies. The mass growth rate is time-dependent, 
slowing at late times (e.g. Wechsler et al. 2002). The cosmo- 
logical context within which clusters evolve may in fact play a 
significant role in their thermal histo ry, as suggeste d by recent 
hydrodynamic sirnulations (e.g.lMot l et al. 2004; Burns et al.l 
l2007tlNagai et al.ll2007tlMcCarthv et al.ii2007h) . 

There are a variety of ways in which the cosmologi- 
cal setting can provide additional heating mechanisms. A 
generic class of such mechanisms may be called "gravi- 
tational" in the sense that the energy available for heat- 
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ing ultimately comes from the gravitational potential en- 
ergy of infalling material. DF heating is one example of 
this class, but there are others. In particular, if small 
gaseous clumps are accreted, they may transfer their grav- 
itational energy to the ICM via ram pressure drag and lo- 
cal shocks (e.g. [Murray & Linll20 04l: iMcCarth v eT5]l2007bt 
iKhochfar &"Ostr iker 2007; Dekel & Birnboim'2007|). Major 
mergers can significantly mix the ICM, effectively transport- 
ing the abundant reservoir of energy in the clus ter outskirts to 
the cooler inner regions (e.g. lBurns et al.ll2007l) . 

Recent observations have suggested that the fraction of 
CC clusters is increasing with time sinc e at least z '--^ 1 
(IVikhUnin et al.ll2006al: lO'Hara et al.ll2007h . If gravitational 
heating plays a significant role in the thermal history of the 
ICM, then one may expect that as the accretion rate slows at 
late times (Wechsleret al. 2002), radiative losses would be- 
come increasingly dominant, in at least qualitative agreement 
with these recent observations. More generally, it is clear that 
the time-evolution of the abundance of observed CC clusters 
will provide strong, unique constraints on the importance of 
various heating mechanisms of the ICM. 

While it may be appealing to invoke major mergers or other 
strong mixing processes as a means to transfer energy into the 
inner cooling regions of clusters, such processes have poten- 
tial draw-backs. The most serious is that observed CC clus- 
ters have strong met allicity gradients within the cooling re- 
gion (iDe Grandi & M olendi 2001). If the metals are produced 
by type la SNe within the central galaxy then the t imescale for 
gener ating the observed gradient is 5 Gyr (Bo hringer et al.l 
|2004'), which is much longer than th e cooling time of the ICM 
in most observed clusters (Sanderso n et al.|20 06). Sedimenta- 
tion may also contribute to the metallicity gradient on relevant 
timescales, but its relevance is much more uncertain because 
it can be highly sup pressed in the presence of even modest 
magnetic fields (e.g. iFabian & Pringlelll977h . If mergers are 
capable of mixing the ICM such that the inner cooling region 
is effectively heated, it seems plausible that they may also de- 
stroy the metallicity gradient. Detailed numerical simulations 
are required to address these issues. 

It is clear that 3D cosmologically embedded simulations 
are preferable to ID simulations of the type presented herein 
because there are potentially relevant physical processes that 
cannot be adequately captured in the latter approach. While 
we have highlighted the potential importance of the cosmo- 
logical context, it is also true that certain physical processes 
cannot be captured simply because our simulations are in ID. 
For example, processes such as convection and the Rayleigh- 
Taylor instability cannot be captured in a ID simulation. Ad- 
ditional insight into the cooling flow problem can thus be 
gained by simply extending the types of simulations we have 
run into 3D. 

In spite of these fact, the motivations for focusing on ID 
simulations are several-fold. First, cosmologically embedded 
simulations are very time-consuming — a parameter space 
study similar to what we have presented in this work is not 
currently possibl e with such simulations. Se cond, DF is a res- 
onant process ( Trem aine & We inberg 1984), and so the parti- 
cle requirement to adequately resolve this phenomenon may 
be considerably higher than what is achievable in the current 
generation of cosmologically embedded hydrodynamic sim- 
ulations. It is not sufficient to resolve well the dark matter 
halo and subhalos. Since much of the DF heating is caused 
by the inspiralling of the much smaller stellar systems, a res- 
olution small compared to the ^10 kpc half-light radii of gi- 



ant ellipticals is required. Thus, while our approach to DF 
heating has been idealized, it was necessary to explore its 
possible effects in this parameterized manner because it is 
not clear that current 3D simulations are adequately resolving 
this phenom enon (see discussion in Faltenba cher et al.ll2005l : 
iNaab et all 12007). Finally, even with sufficient resolution, 
there are a number of physical effects that we have demon- 
strated are potentially relevant for the thermodynamics of the 
ICM that are not included in the majority of cosmologically 
embedded hydrodynamic simulations. Such phenomena in- 
clude thermal conduction, which ultimately requires the in- 
clusion of magnetic fields, type la supernovae energy injec- 
tion, which is included in some but not all current simulations, 
and energy injection from accretion onto a central black-hole, 
which we now discuss in detail. 

4.2. A Self-Regulating Heating Mechanism 

In the present work we have focused on heating processes 
that are not generally thought of as "feedback" mechanisms. 
While DF heating is in principal a feedback mechanism since 
the heating becomes inefficient at very and low high Mach 
numbers, in practice this feedback is too weak to prevent run- 
away heating or cooling. The preceding discussion of the 
importance of the cosmological context suggests that addi- 
tional gravitational energy is available that can contribute to 
the heating of the ICM. However, the incredibly short cool- 
ing times of many CC clusters (<C 1 Gyr) indicates that heat- 
ing from without, via gravitational processes, must still be 
fine-tuned in order to prevent runaway cooling. Heating from 
within, via the release of energy from accretion onto a black 
hole, is in many ways more appealing because such a process 
is explicitly self-regulating. 

Black hole accretion-mediated feedback mechanisms 
(e.glC iotti & Ostriker 2001; Ruszkowski & Begelman 200l 
'Kaiser & Binney, i2003c iGuo & Ohi i2007i; iCiotti & Ostrikerl 
2007) are qualitatively different from the mechanisms we 
have explored. These self-regulating feedback processes are 
thought to work schematically as follows. As the gas begins 
to cool in the inner regions, its thermal pressure no longer 
provides sufficient support to the overlying material, and thus 
matter flows inward and eventually accretes onto the central 
black hole. The mechanical energy provided by AGN activ- 
ity, which is proportional to the rate of mass infall, then heats 
up the surrounding gas until the thermal pressure, which is 
now increasing due to the energy injection, halts the inward 
flow of matter, thereby diminishing the accretion-driven en- 
ergy injection. The cycle then begins again in a regulatory 
fashion, and the ICM thus neither heats nor cools catastroph- 
ically. It is this regulatory feedback processes which seems 
most promising in explaining the properties of observed intr- 
acluster media because the relevant physical parameters need 
not be tuned to any particu lar values for hea ting to balance 
cooling generically (see e.g. lGuo & Ohll2007l) . 

A serious constraint on AGN-related feedback is that it may 
destroy the observed metalUcity gradient in CC clusters if the 
energy deposition at t he cluster center i s suffi cient to drive 
convection. Recently, IVoit & Donahud (l2005h have shown 
that AGN feedback can be effective at balancing radiative 
cooling and maintaining the observed metallicity gradient if 
the AGN outbursts are rather gentle, occur every ~ 10^ years 
and last in duration for 10^ years. It remains to be seen if 
this proposal can be confirmed with both numerical simula- 
tions and, ultimately, observations. 

While AGN heating may be a promising candidate at gen- 
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erating and maintaining intracluster media in thermal bal- 
ance over cosmological time-scales, simulations of the ICM 
which include it do so only in a rather simplistic fashion, 
in part because a detailed understanding of the physical pro- 
cesses involved in black hole accretion currently eludes us. 
The way in which energy is transfered from the black hole 
to the surrounding ICM is also obscure. Until these and 
related issues are better understood, we should continue to 
seek out other potential self-regulating heating processes that 
could potentially be relevant for the thermodynamics of the 
ICM. Nonetheless, varied obser vations of sound waves, bub- 
bles, a nd strong radio ern ission (iBirzan et al.ll2004HBest et all 
l2005t iFabian et al.ll2006l) . predominantly in CC clusters, in 
additi on to their very short cooling times feanderson et al.l 
12006'), strongly suggests that AGN-related activity plays an 
important role in the thermodynamics of the ICM. 

5. SUMMARY 

We have presented the results from a series of ID simula- 
tions aimed at understanding the importance of type la super- 
novae heating, thermal conduction and DF heating (due to the 
orbital motions of satellite galaxies) on the thermal proper- 
ties of the intracluster media of clusters. Both initially non- 
cooling core (NCC; i.e. isothermal) and cooling core (CC) 
clusters were simulated, for a wide range of initial central 
electron densities and gas temperatures. 

Marginalizing over the uncertain efficiencies of DF heating 
and thermal conduction, it is clear that only NCC clusters with 
central electron densities nj,,o ^ 0.02 cm~^ can be maintained 
in thermal balance over a Hubble time if both DF and conduc- 
tion operate; neither mechanism alone can generate generi- 
cally stable ICM at these densities. At higher densities no 
reasonable amount of conduction or DF heating can prevent 
runaway cooling in NCC clusters. 

The failure of the combination of conduction and DF heat- 



ing at either generating or maintaining observed intracluster 
media is more pronounced for CC clusters, which have ob- 
served temperature profiles that decline by a factor of ~ 2 - 3 
from the outer to inner regions. This temperature drop is ex- 
tremely difficult to maintain in the face of thermal conduction 
in our simulations, because thermal conduction acts to erase 
temperature gradients. In fact, of the 0(10^) simulations run 
with a wide range of DF heating and conduction efficiencies 
and initial gas densities, none generated stable CC cluster pro- 
files. Sin ce CC clusters co nstitute ^ 70% of observed clus- 
ters (e.g. iPeres et al.|[T998h . we regard this failure as strong 
evidence that other heating processes in the ICM must be at 
work besides conduction and DF heating. 

Our results demonstrate that there are numerous energy 
reservoirs capable of supplying enough energy to offset ra- 
diative losses. The crux of the cooling flow problem therefore 
lies not in finding one or more mechanisms capable of pro- 
viding enough energy to the ICM, but rather in finding one 
or more mechanisms that can supply the energy in a way that 
maintains thermal equilibrium in the ICM over cosmologi- 
cal time-scales. Only low-density NCC clusters can be main- 
tained with DF heating and conduction. The processes ex- 
plored herein are not manifestly self-regulating, and they thus 
require fine-tuning in order to generate the observed proper- 
ties of CC and high-density NCC intracluster media generi- 
cally. For these types of observed intracluster media, it seems 
likely that an explicitly self-regulatory feedback process such 
as black hole accretion-powered energy injection (i.e. AGN 
feedback) is required. 
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APPENDIX 

A. THERMAL STABILITY FOR GENERIC HEATING MECHANISMS IN THE PRESENCE OF A RELATIVISTIC FLUID 

The generalized Field criterion (lFieldlll965t lBalbus|[T986l) states that a gas is thermally unstable to isobaric perturbations if: 

d{C/T) 



dT 



<0, 



(Al) 



where £ is the net loss function per unit mass defined such that pgC = A - F where A and F are the cooling and heating rates per 
unit volume. In what follows the subscript "g", denoting gaseous quantities, will be omitted for brevity. 
The following identity holds for a general gas in thermal equilibrium {pC = 0): 
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P dp 
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The first and third derivatives on the right-hand side are determined by the relevant heating and cooUng mechanisms while the 
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second is set by the various sources of pressure support. The following identity will be useful: 
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Assume that the pressure support is provided by both the thermal pressure of the gas and a relativistic fluid: 

P = Px^+Pre\ = KipT+P,,x(p), 

where K\ is a constant and Piei is a function only of density. It then follows that 

dP 
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and 
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and hence: 
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where in the last equality we have defined a' = Pth/Aei- 

Now assume that cooling is dominated by thermal Bremsstrahlung radiation and that the heating has a simply power-law 
dependence on p and T. Then: 

p/: = Aop^r'/2-rop'3r-'/2+'5, (as) 

where Aq and Fq are constants and /3 and S are the power-law indices for heating, defined in a way that will prove useful below. 
After some algebra, and using the fact that pC = 0, we have 



d(p£) 
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Combining these with Equations I A2 1 and I A7 1 we have, after more algebra 
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The gas is thus thermally unstable if: 

6> 



d\np \t 



a' + 7 



l-a I 
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(A12) 



where we have assumed that Prei « p'^ and d efined a = Pjei/A ot- 

This equation recovers the result found in lKim et al.l (l2005h . There they explored the stability of DF h eating, where P = 1, in 
the absence of relativistic pressure (a = 0). The instability criterion in this case is 5 > 0. As shown in iKim et al.l (l2005h . DF 
heating generally yields a temperature dependence such that 6 >0, and thus DF heating alone is thermally unstable. 

This equation also implies that heating sources which scale as p^ are thermally stable if the temperature dependence scales as 
a power less than 1 /2, which is physically intuitive since p^T^I^ is the scaling of the cooling function. Moreover, if the heating 
source scales as then the relativistic fluid has no influence on the thermal stability. 



